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Abstract 



By comparing structure functions measured in neutrino and charged lepton 
deep inelastic scattering, one can test the validity of parton charge symmetry. 
New experiments allow us to make such tests, which set rather tight upper 
limits on parton charge symmetry violation [CSV] for intermediate Bjorken 
x, but which appear to show sizable CSV effects at small x. We show that 
neither nuclear shadowing nor contributions from strange and antistrange 
quark distributions can account for the experimentally observed difference 
between the two structure functions. We are therefore forced to consider the 
possibility of a large CSV effect in the nucleon sea quark distributions. We 
discuss the consequences of this effect for other observables, and we propose 
an experiment which could detect a large CSV component in the nucleon sea. 







I. INTRODUCTION 



In discussing the strong interaction, it is customary to assume the validity of charge 
symmetry, which interchanges protons and neutrons (simultaneously interchanging up and 
down quarks). For example, all phenomenological analyses of deep inelastic scattering data 
in terms of parton distribution functions assume charge symmetry from the beginning. Our 
faith in charge symmetry is justified from our experience in nuclear physics, where this 
symmetry is respected to a high degree of precision. Most experimental low-energy tests 
of charge symmetry find that it is good to at least 1% in reaction amplitudes Until 
recently such an assumption seemed to be justified, as there was no compelling experimen- 
tal evidence against parton charge symmetry. The quantitative evidence which could be 
extracted from high energy experiments, although not particularly precise, was consistent 
with charge symmetric parton distributions 0. 

Experimental verification of charge symmetry is difficult, partly because the relative 
charge symmetry violation (CSV) effects are expected to be small, requiring high precision 
experiments to measure CSV effects, and partly because CSV often mixes with parton flavor 
symmetry violation (FSV). Recent experimental measurements by the NMC Collaboration 

demonstrating the violation of the Gottfried sum rule ||, have been widely interpreted 
as evidence for what is termed SU(2) FSV. The measurement of the ratio of Drell-Yan cross 
sections in proton-deuteron and proton-proton scattering, first by the NA51-Collaboration at 
CERN H and more recently by the E866 experiment at FNAL |7J , also indicate substantial 
FSV. However, both of these experiments could in principle be explained by sufficiently large 
CSV effects HH, even in the limit of exact flavor symmetry. In view of these ambiguities 
in the interpretation of current experimental data, it would be highly desirable to have 
experiments which separate CSV from FSV. A few experiments have been already proposed 
|iO| , [TTH and could be carried out in the near future. 



Recent experiments now allow us for the first time to make precision tests which could 
put tight upper limits on parton CSV contributions. The NMC measurements of muon DIS 
on deuterium [12\ provide values for the charged lepton structure function F^x^Q 2 ). In a 
similar Q 2 regime the CCFR Collaboration [IB| extract the structure functions F£(x, Q 2 ) 
from neutrino-induced charge-changing reactions. As we show in sec. II, the "charge ratio", 
which can be constructed from these two quantities (plus information about the strange 
quark distribution) can in principle place strong constraints on parton CSV distributions. 
We will show that, for intermediate values 0.1 < x < 0.4, the agreement between the two 
structure functions is impressive, and provides the best upper limit to date on parton CSV 
terms. However, the charge ratio shows a substantial deviation from unity in the region 
x < 0.1, which might suggest surprisingly large charge symmetry violation. In a recent 
Letter [0] we argued that the data supported this conclusion. However, several important 
corrections have to be applied to the data before any conclusions can be reached. These 
corrections are especially important for the neutrino cross sections. 

In sec. II we discuss the uncertainties involved in the analysis of the data. Most correc- 
tions have already been accounted for in the present experimental analysis. We particularly 
focus on two aspects of the neutrino reactions: heavy target corrections and effects due 
to strange and antistrange quark distributions. In sec. Ill we demonstrate that neither of 
these effects are sufficient to account for the apparent discrepancy at small x. The charge 
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symmetry violating distributions can be obtained from a combination of neutrino charged 
current structure functions, muon structure functions and strange quark distributions ex- 
tracted from dimuon production in neutrino reactions. We construct such a combination 
and extract the CSV terms. Assuming the validity of the experimental data, we find CSV 
effects on the order of 25% of the sea quark distributions at low x. In sec. IV we discuss the 
consequences of such large CSV effects on other observables. We examine the role played 
by CSV in the extraction of the FSV ratio d/u, in the Gottfried sum rule and in the exper- 
imental determination of the Weinberg angle sin 2 8w In sec. V we suggest an experiment 
which could measure the substantial CSV suggested by our analysis. 

II. COMPARING STRUCTURE FUNCTIONS FROM NEUTRINO AND 
CHARGED LEPTON REACTIONS 

Our analysis of parton charge symmetry violation is based on the "charge ratio," which 
we review here. This depends on the ratio of F2 structure functions extracted from charged 
lepton reactions with those from neutrino charge-changing reactions. Because neutrino 
cross sections are so small, at present the structure functions can only be measured for 
heavy targets such as iron. Furthermore, in order to obtain useful statistics, the data must 
be integrated over all energies for a given x and Q 2 bin. As a result, only certain linear 
combinations of neutrino and antineutrino structure functions can be obtained. The process 
by which we attempt to extract parton CSV contributions is complicated, and requires input 
from several experiments. In this section we review this process in detail. 



Structure functions measured in neutrino and muon deep inelastic scattering are inter- 
preted in terms of parton distribution functions. Since the operation of charge symmetry 
maps up quarks to down quarks, and protons to neutrons, at the level of parton distribu- 
tions, charge symmetry implies the equivalence between up (down) quark distributions in 
the proton and down (up) quark distributions in the neutron. In order to take CSV in the 
parton distributions into account, we define the charge symmetry violating distributions as 



where the superscripts p and n refer to quark distributions in the proton and neutron, 
respectively. The relations for CSV in antiquark distributions are analogous. If charge 
symmetry were exact then the quantities 5u(x) and 8d(x) would vanish. 

In the quark-parton model the structure functions measured in neutrino, antineutrino and 
charged lepton DIS on an isoscalar target, N , are given in terms of the parton distribution 
functions and the charge symmetry violating distributions defined above by 



A. The "Charge Ratio" and Charge Symmetry Violation 



Su(x) = u p (x) - d n (x) 




(1) 



F% N °(x,Q 2 ) = x[u{x) + u{x) + d(x) + d{x) + 2s(x) + 2c{x) - 5u{x) - 8d(x)\ 
F* No (x, Q 2 ) = x[u{x) + u{x) + d(x) + d(x) + 2s{x) + 2c(x) - 5d{x) - 5u(x)} 
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xF^ No {x, Q 2 ) = x[u(x) + d(x) - u(x) - d(x) + 2s(x) - 2c(x) - 5u(x) + 6d(x)] 
xF* N °(x, Q 2 ) = x[u(x) + d(x) - u(x) - d(x) - 2s(x) + 2c(x) - 5d(x) + 5u(x)\ 

5 2 8 

F% N °(x,Q 2 ) = — x[u{x) +u{x) +d(x) + d(x) + -(s(x) + s(x)) + -(c(x) + c(x)) 

4 1 

- -(5d(x) + Sd{x)) - ~{Su{x) + Su(x))] (2) 

D 

Here, and in the following, quark distributions without superscripts denote quark distribu- 
tions in the proton. From now on, we will disregard charm quark contributions to the struc- 
ture functions. Since phenomenological parton distribution functions assume the validity of 
charge symmetry, possible CSV effects are folded into the commonly used phenomenologi- 
cal parton distribution functions in a highly non-trivial way. Nevertheless, using the above 
relations, it is possible to test the validity of charge symmetry by building appropriate lin- 
ear combinations or ratios of the measured structure functions. One such possibility is to 
calculate the "charge ratio", which relates the neutrino structure function to the structure 
function measured in charged lepton deep-inelastic scattering 



|F*(.,«'|-i(»(') + i("))/f 

s(x) — s(x) 45u(x) — 8u(x) — 45d(x) + Sd(x) 



Q(x) 5Q(x) 



(3) 



Here, we defined Q(x) = J2 q =u,d,s{Q{ x ) + < 7( :r )) — 3(s(sc) + s(x))/5, and we have expanded Eq. 
|3] to lowest order in small quantities. From Eq. |3] we see that any deviation of the charge 
ratio from unity, at any value of x, would be due either to CSV effects or to different strange 
and antistrange quark distributions. Analogous relations could be obtained using structure 
functions from antineutrinos, or from a linear combination of neutrino and antineutrino 
structure functions. For example, we can derive 



n c (x,Q 2 ) EE 



F£ No {x,Q'' 



F% No (x,Q 2 ) - x(s(x) + s(x))/6 
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3 (5u(x) + Su(x) — 5d(x) — 8d(x) ) 
«H ^ = '-. (4) 

In Eq. | Ti N °{x, Q 2 ) = (FZ No (x, Q 2 ) + F^ N °(x, Q 2 ))/2 is the average of the structure func- 
tions from neutrino and antineutrino reactions; deviations from one in the ratio 1Z c (x) depend 
only on parton CSV contributions, and have no contribution from strange or antistrange 
quark distributions. 



The recent measurement of the structure function i 7 ^ by the CCFR-Collaboration fl3 
makes it possible to carry out a precise comparison between F%{x, Q 2 ) and F$(x, Q 2 ) for the 
first time. The CCFR-Collaboration compared the neutrino structure function F£(x, Q 2 ) 
extracted from their data on an iron target |13|] with F^x, Q 2 ) measured for the deuteron 



by the NMC Collaboration In the region of intermediate values of Bjorken x (0.1 < x < 
0.4), they found very good agreement between the two structure functions. In this x region, 
this allows us to set upper limits of a few percent on parton CSV contributions. On the 
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other hand, in the small rc-region (x < 0.1), the CCFR group found that the two structure 
functions differ by as much as 10-15%. This can be seen in Fig.|l| where the "charge ratio" has 
been obtained by integrating over the region of overlap in Q 2 of the two experiments. The 
open and solid circles in Fig. |l] represent two different ways of calculating nuclear shadowing 
corrections, as we will discuss later. 



B. Extracting Structure Functions From Neutrino Cross Sections 

In order to perform tests of parton distributions through, say, the charge ratio of Eq. |3|, 
we need the structure functions from neutrino charge-changing reactions on a free proton 
and neutron. These are written in terms of parton distributions in Eq. ^. Because of 
the extremely small cross sections for neutrino-induced reactions, we are able to obtain 
statistically meaningful cross sections only from heavier targets such as iron. We then 
have to make the following corrections in order to extract the neutrino structure functions 
on "free" nucleons, averaged over proton and neutron, Fp uN °(x, Q 2 ), (j = 2,3): i) The 
nuclear structure functions Fp uFe (x,Q 2 ) must be extracted from the cross sections; ii) The 
nuclear structure functions need to be corrected for the excess of neutrons in iron (isoscalar 
effects); iii) Kinematic corrections must be applied to account for heavy quark thresholds, 
particularly charm quark threshold effects (Eq. |3| is valid only well above all heavy quark 
thresholds); iv) Heavy target corrections must be applied, to convert structure functions 
for nuclei to those for free protons and neutrons; v) The neutrino and muon cross sections 
must be properly normalized. In order to test charge symmetry, all these corrections have 
to be taken into account. The data have already been corrected for normalization, isoscalar 
and charm threshold effects by the CCFR-Collaboration in their analyses There is a 



thorough discussion of these points in the thesis by W. Seligman |Tj|. Here we will review 
how the nuclear structure functions are extracted from the cross sections, the heavy target 
corrections for neutrino reactions, and the role of both strange quarks and CSV effects in 
neutrino structure functions. 

The cross sections for neutrino and antineutrino scattering on a nuclear target containing 
A nucleons can be written as 

IxlQ 2 = ^2 {F * PA{X > Q2)±XF3U ' PA{X > Q2)) + j( F 2 PA ^Q 2 )TxFr A (x,Q 2 ))]. (5) 

In Eq. |5] the upper (lower) sign is associated with neutrino (antineutrino) cross sections. 
We have assumed the validity of the Callan-Gross relation and neglected terms of order 
Mxy/2E, and we introduced the variable £ = (1 — y). It would be straightforward to remove 
these assumptions. With a large enough count rate, the x and y dependence of the cross 
sections could be separately measured. By plotting the measured differential cross sections 
for fixed x and Q 2 as a function of £ 2 , the structure functions F2 and F3 can be determined 
from the slopes and intercepts of the resulting straight lines. The crucial question is, of 
course, whether the statistics of the experiment are sufficient for the structure functions to 
be extracted in this way. 

To illustrate this problem we calculated the statistical errors in each energy bin. For 
this calculation, we used the experimental determined fluxes, the total and differential neu- 
trino and antineutrino cross sections to obtain the expected number of events in a given 
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x, Q 2 and energy bin. We estimated the statistical errors using Act = o/yN. In Fig|J 
<t u,u (x,Q 2 ,!; 2 )/(G 2 7 /2tix) is plotted as a function of £ 2 . The solid lines are the results us- 
ing the CTEQ parton distribution functions and assuming the validity of the Callan-Gross 
relation. The dotted lines are the results obtained without using the Callan-Gross rela- 



tion. Here, we used the parametrization of Whitlow [|T(J for the ratio of the longitudinal 
and transverse photo-absorption cross sections. The current statistics do not allow one to 
extract the individual structure functions. The error bars represent the expected statistical 
errors. An order of magnitude more events would be necessary to decrease the statistical 
errors sufficiently that one could consider extracting the structure functions directly, and 
systematic errors would further complicate this analysis. 

Since the number of events is so small that individual structure functions cannot be 
extracted from the data, the cross sections in a given x and Q 2 bin are integrated over all 
energies. After this integration is performed, Eq. [| can be written as two linear equations, 
one for neutrino and the other for antineutrino events: 

N v (x, Q 2 ) = A£ Ff e (x, Q 2 ) + A3 xFf e (x, Q 2 ) 

N*(x,Q 2 ) = A*F* Fe (x,Q 2 )-A*xFf e (x,Q 2 ) . (6) 

In Eq. |6| N u (N u ) is the number of neutrino (antineutrino) events in a given x and Q 2 - 
bin integrated over the incident neutrino and antineutrino energies. A\ and A\ (i = 2, 3) 
represent the coefficients, Ai(y), of the structure functions multiplied by the neutrino and 
antineutrino fluxes, <& V (E) and <$> U (E), respectively, and integrated over all energies 

A» = J dEAi(y) &{E) 

A\= [ dEA(y)^(E). (7) 



The individual structure functions for neutrino and antineutrino reactions are extracted 
by taking linear combinations of the relations in Eq. |] and making corrections using phe- 
nomenological parton distribution functions. For example, from Eq. we see that for an 
isoscalar target, F^ (x,Q 2 ) = F% °(x, Q 2 ) if charge symmetry is valid and s(x) = s(x). 
Thus we can form linear combinations of the terms in Eq. ^| such that these terms cancel 
and we are left only with the F3 structure functions. Similarly, assuming charge symmetry 
we have 

F» N °{x,Q 2 ) - F£ N °(x,Q 2 ) = 2[s(x) - s(x)) . 

We can then take a linear combination of the terms in Eq. |] which gives this function. If 
the strange quark distribution is taken from a phenomenological model, we can extract a 
linear combination of the F 2 structure functions for neutrinos and antineutrinos on a nuclear 
target. 

We will discuss how the structure functions are extracted, and particularly the role of 
CSV and strange quark distributions in this process. However, at this stage we review how 
heavy target corrections are calculated, in order to extract the structure functions for free 
nucleons from those measured on a heavy nuclear target. 
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C. Heavy Target Corrections in Neutrino Reactions 



As is well known, the structure functions measured on heavy targets are not equal to those 
observed for light targets such as the deuteron. At small x values, nuclear shadowing effects 
play a major role; at large x, nuclear Fermi motion effects dominate, and at intermediate x 



"EMC" effects play a significant role fl7| . Such effects have been systematically measured 
in charged lepton reactions. 

In analyzing neutrino scattering data, it is generally assumed that heavy target correc- 
tions will be the same as those observed in charged lepton reactions. A priori, there is no rea- 
son to assume that neutrino and charged lepton heavy target corrections should be identical. 
Heavy target corrections for neutrinos are generally applied by multiplying the experimental 
structure functions at a given x value by the quantity R = F2 (x, Q 2 )/F| r> (a;, Q 2 ), the ratio 
between the F2 structure function measured on heavy targets and that of the deuteron for 
charged lepton deep inelastic scattering, at the same x value. However, as is well known, 
shadowing corrections are very much Q 2 dependent for smaller Q 2 values (where a consid- 
erable part of the available data was taken), and the Q 2 and x-dependence of the data are 
strongly correlated because of the fixed target nature of these experiments. 

We re-examined heavy target corrections to deep-inelastic neutrino scattering, focusing 
on the differences between neutrino and charged lepton scattering and on effects due to the 
Q 2 -dependence of shadowing for moderately large Q 2 . This work will be published elsewhere 
]TE| ; here we briefly review the results of that work. We used a two phase model which has 
been successfully applied to the description of shadowing in charged lepton D1S p!9| , |20| . 
In this approach, vector meson dominance is used to describe the low Q 2 virtual photon 
or W interactions, and Pomeron exchange is used for the approximate scaling region. In 
generalizing this approach to weak currents, the essential differences in shadowing between 
neutrino and charged lepton deep inelastic scattering are: (i) the axial-vector current is only 
partially conserved, in contrast to the vector current; and (ii) the weak current couples not 
only to vector but also to axial vector mesons f21|-|24||. 

Partial conservation of the axial current (PCAC) requires that the divergence of the axial 
current does not vanish but is proportional to the pion field for Q 2 = 0. This is Adler's 
theorem [p5|1 , which relates the neutrino cross section to the pion cross section on the same 
target for Q 2 = 0. Thus, for low Q 2 (~ m 2 ) shadowing in neutrino scattering is determined 
by the absorption of pions on the target. For larger Q 2 - values the contributions of vector and 
axial vector mesons become important. The coupling of the weak current to the vector and 
axial vector mesons and that of the electro-magnetic current to vector mesons are related 
to each other by the 'Weinberg sum rule" f 2 + = f 2 ^ = 2f 2 . Since the coupling of the 
vector (axial vector) mesons to the weak current is twice as large as the coupling to the 
electro-magnetic current, but the structure function is larger by a factor of ~ 18/5 in the 
neutrino case, we expect that shadowing due to VMD in neutrino reactions is roughly half 
of that in charged lepton scattering. 

For larger Q 2 -values, shadowing due to Pomeron exchange between the projectile and two 
or more constituent nucleons dominates. Since Pomeron-exchange models the interaction 
between partons in different nucleons and the scattering of the W takes place on only one 
parton, this processes is of leading twist in contrast to the VMD and pion contributions. 
The coupling is given by the coupling of the photon or W to the quarks in the exchanged 
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Pomeron. It changes in the same way as the structure function does in switching from 
neutrino to charged lepton scattering. Thus, for large Q 2 values (> 10 GeV 2 ), shadowing in 
both cases should have approximately the same magnitude. In the intermediate Q 2 -region 
(1 < Q 2 < 10 GeV 2 ), where VMD is relatively important, we expect to see differences 
between shadowing in neutrino and charged lepton scattering. We recall that this is precisely 
the region where the discrepancy between CCFR and NMC is significant. There are also 
nuclear effects in the deuteron. However, because of the low density of the deuteron, these 
are (relatively speaking) very small and have a negligible effect on the charge ratio. 

We calculated the shadowing corrections to the CCFR neutrino data using the two-phase 
model of Ref . p9|j20]] . With this corrected CCFR data, we calculated the charge ratio R c of 



Eq. [| between CCFR and NMC data. The result is shown in Fig.|l|. The open triangles show 
the charge ratio when no shadowing corrections are used. The open circles show the charge 
ratio when heavy target shadowing corrections from charged lepton reactions are applied 
to the neutrino data, and the solid circles show the result when the neutrino shadowing 
corrections from our two-phase model are applied. At small x, using the "correct" neutrino 
shadowing corrections reduces the deviation of the charge ratio from unity. Nevertheless, 
the charge ratio is still not compatible with one at small x. In summary, properly accounting 
for shadowing corrections in the neutrino structure function decreases, but does not resolve, 
the low- a; discrepancy between the CCFR and the NMC data. 



D. Strange Quark and CSV Contributions to Structure Functions 

In Eq. ^| we showed that, after integrating neutrino charged-current cross sections over 
all energies, we obtain two equations in four unknowns, the structure functions F 2 and F 3 
for neutrino and antineutrino reactions. If the neutrino and antineutrino structure functions 
were equal, F% Fe (x,Q 2 ) = Fg Fe (x,Q 2 ), with an analogous relation for xFg Fe (x,Q 2 ), then 
Eq. H would provide two linear equations in two unknowns. As we discussed previously, 
several corrections need to be applied before we can extract the structure functions on a 
"free" isoscalar target N , and compare the structure functions to the parton distributions 
given in Eq. £| First, since iron is not an isoscalar target we need to make corrections 
for the excess neutrons. Next, we need to estimate the contributions from strange quark 
distributions and charge symmetry violating parton distributions. Finally, we need to make 
heavy target corrections as reviewed in the preceding section. 

We begin by splitting the neutrino and antineutrino structure functions on iron into 
isoscalar and non-isoscalar parts. For a target with Z protons and N = A — Z neutrons we 
define the quantity (3 = (N — Z)/A: 

pu,DFe _ }^^pv,up _|_ pu,un-^ _ @_^pv,up _ pu,Dn^ 

The first term on the right of Eq. |8] corresponds to the neutrino and antineutrino structure 
functions on an isoscalar target, N . The second terms include corrections arising from the 
non-isoscalarity of the target. In the absence of CSV, these corrections are basically given 
by the difference between up and down valence quark distributions and have been taken into 
account in the extraction of the structure functions. However, the non-isoscalarity of the 
target leads also to CSV corrections. 
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We define the sum and difference of the neutrino and antineutrino structure functions 
on a target A as 

ATt = ±[F? A -FT A ) ; (9) 
the structure functions p^ vFe can then be written as 

F u,PFe = ^N, ± AJ7 N -^{[ff- ± [Aj 7 ? - A J 7 ™] } . (10) 

Here, "+" and "— " refer to the neutrino and antineutrino structure functions, respectively. 
The last three terms of the right hand side of Eq. [10] contain corrections coming from excess 
neutrons, strange quarks, CSV and s(x) ^ s(x). Correcting the data for excess neutrons 
and for strange quark contributions corresponds to subtracting the number of events due to 
the corresponding corrections from the left hand side of Eqs. |6] 



i=2 



i=2 



In Eq. [0], we have calculated corrections to the structure functions from excess neutrons 
and strange quarks, and have used these to produce the effective number of events on 
the left hand side of Eq. [11]. (&Fj) n|S and (SFAcsv Te ^ eT t° corrections arising from excess 
neutrons, strange quark distributions, because of charge symmetry violation and s(x) ^ s(x), 
respectively. The CCFR Collaboration assumed the validity of charge symmetry, and they 
also took s(x) = s(x) based on the results of a next to leading order [NLO] analysis of 



dimuon production in neutrino-induced reactions [26]. We have left the correction terms 
coming from CSV and s(x) ^ s(x) on the right hand side of Eq. |6| as these have been 
absorbed into the extracted structure functions. Under the assumption of charge symmetry 
and s(x) = s(x), Eq. [11] simplifies to 



N v - £ AttSFZUs = A" jtCcfr,a + A u xT ccfra 

i=2 

3 

N p - X)(-1)M? {5F») n , s = A\ j^ CFR > A - A\ xT^ CFRA . (12) 

8=2 

These equations provide a system of two linear equations for the two nuclear structure 
functions j?ccfr,a an( ^ j^ccfr^a^ p rom these structure functions we can calculate the 

structure functions for a "free" nucleon target using the heavy target correction factors 
described in the previous section. The resulting structure functions Tf CFR still contain 
charge symmetry violating contributions and terms proportional to s(x) — s(x), as can be 
seen from Eq. ||. To relate the measured structure functions, Ff CFR to the various parton 
distributions, we take the sum and difference of the measured number densities in Eqs. O 
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and [12] (for a fixed energy) and compare the coefficients of Ai(y). In this way, we see that 
the measured structure functions, Ff , can effectively be identified with a flux weighted 
average of the neutrino and antineutrino structure functions, F^ N ° and F" N ° } and correction 
terms arising from CSV effects 

pCCFR = T N + (2a _ 1)AJ7 N -fL[J*- T-\CSV 

- {2a ~ 1)f3 [AF?-AF?] C sv. (13) 

Here, we defined the relative neutrino flux, a, as a = $> v /(<& v + $ I/ ). The experimental value 
of a depends on the incident neutrino and antineutrino energies and is also different for the 
E744 and E770 experiments. Because of the kinematical constraint y < 1, relative fluxes at 
energies > 150 GeV are relevant for small x. Here, a ~ 0.83 |13[] so that F 2 FR (x, Q 2 ) can 
be approximately regarded as a neutrino structure function. 

The different contributions to F2 can be expressed in terms of the quark distribution 
functions 

\\n - n\csv = -\^2°}csv = I [6u(x) + 5u(x) + 5d{x) + 5d{x)) 
i[A^f - AF2} CSV = -| [Sd(x) - 5d(x) - 5u(a;) + 5u(x)} 

A jrN = -l[AJ* - A^] CSV + x[s(x) - s(x)\. (14) 

The second expression in Eq. |14| is obtained by subtracting the F 2 structure function for 
neutrinos on protons from that for antineutrinos on protons; from this is subtracted the 
corresponding term for neutrons. It depends only on charge symmetry violation in the 
valence quark distributions. The last expression in Eq. [14| is obtained by taking the difference 
between neutrino and antineutrino F2 structure functions on an isoscalar system. It also 
depends on valence quark CSV, and has an additional contribution from the difference 



between strange and antistrange parton distributions. The first term in Eq. TA is obtained 
by averaging the F2 structure functions over neutrino and antineutrino reactions, and taking 
the difference of the F2 structure functions measured on proton and neutron targets. This 
quantity is free from strange quark effects, and is also sensitive to CSV in the sea-quark 
distributions. 



III. EVIDENCE FOR LARGE CHARGE SYMMETRY VIOLATION IN PARTON 

SEA QUARK DISTRIBUTIONS 

The most likely explanation for the discrepancy in the small-x region of the charge ratio 



involves either differences between the strange and antistrange quark distributions f27H30 
or charge symmetry violation. First, we will examine the role played by the strange and 
antistrange quark distributions. Assuming that charge symmetry is exact, the strange and 
antistrange quark distributions are given by a linear combination of the structure functions 
measured in neutrino and in muon DIS, as can be seen from Eqs. 13, O and [14], 
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\f° cfr {x, Q 2 ) - 3F 2 NMC (x, Q 2 ) = \x [s{x) + s{x)\ + §(2a - 1) x [s{x) - s{x)\. (15) 
o 2 o 

Under the assumption that s(x) = s(x), this relation could be used to extract the strange 
quark distribution. However, as is well known, the strange quark distribution obtained in 
this way is inconsistent with the distribution extracted from independent experiments. 



A. Direct Measurement of Sea Quark Distributions 

The strange quark distribution can be determined directly from opposite sign dimuon 
production in deep inelastic neutrino and antineutrino scattering. To leading order in a 
charge- changing reaction, the incoming neutrino (antineutrino) emits a muon and a virtual 
W boson, which scatters on an s or d (s or d) quark, producing a charm (anticharm) quark 
which fragments into a charmed hadron. The semi-leptonic decay of the charmed hadron 



produces an opposite sign muon. The CCFR Collaboration performed a LO j31| and NLO 
analysis |26| of their dimuon data using the neutrino (antineutrino) events to extract the 
strange (antistrange) quark distributions. Their result differs substantially from the strange 
quark distribution extracted from Eq . (|15D , as mentioned above. 

In the dimuon data one extracts the strange and antistrange quark distributions from 
the neutrino and antineutrino data separately. The analysis performed by the CCFR Col- 
laboration suggests that, while there is a difference between the strange and antistrange 
distributions in LO analysis |]3T| they are equal within experimental errors in NLO p6 



However, since the number of antineutrino events is much smaller than that of the neutrino 
events, the errors of this analysis are inevitably large. 

Since the dimuon experiments are carried out on an iron target, shadowing corrections 
could modify the extracted strange quark distribution, and might account for some of the 
discrepancies between the two different determinations of the strange quark distributions. 
The CCFR-Collaboration normalized the dimuon cross section to the "single muon" cross 
section and argued that the heavy target correction should cancel in the ratio. However, 
the charm producing part of the structure function F^^x, Q 2 ) could be shadowed differently 
from the non-charm producing part i^ cp (x, Q 2 ), unless charm threshold effects cancel in 
the shadowing ratio. This could be the case, because vector mesons with higher masses are 
involved in the charm producing part, and because charm production threshold effects have 
to be taken into account in the Pomeron component as well. 

We calculated the shadowing ratio, R = F£ (x, Q 2 ) / F^i^x, Q 2 ), between the structure 
functions on a heavy target and on a deuteron target for both the charm and non-charm 
producing part of the structure function. We took charm production threshold effects into 
account in the Pomeron component through the slow rescaling mechanism by replacing xip, 
which is the momentum fraction of the Pomeron carried by the struck quark, by £/p = 

2 

xjp{l + ). Here m c is the mass of the charm quark. In the VMD component of F^ix, Q 2 ) 
we included the vector mesons D* + (2010), _D* + (2110) and the axial vector partner .0^(2535) 



of D* + [32], which describe the lightest coherent states of the cd and cs fluctuations of the 
H /+ -boson. They have the same coupling to W + as p + and al but have much heavier masses. 
(The cd fluctuations are suppressed by sin 2 O c .) Because of the larger mass of the charmed 
vector mesons (~ 2.5 GeV), we applied a cut at M\ > 6.3 GeV 2 in the diffractively produced 
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invariant mass of the Pomeron component. This is to be compared with M\ > 1.5 GeV 2 in 
the non-charm producing part of the structure function (these cuts are necessary to avoid 
double counting.) Because of the light quark component of the D- mesons, we expect that the 
D-meson-nucleon total cross sections are comparable to the corresponding cross sections of 
lighter mesons with the same light quark content. We use ctd*n ~ ct p n and (Jd*n ~ &<t>N- The 
calculated ratios, R = F 2 A (x,Q 2 )/F 2 D (x,Q 2 ), are shown in Fig. § for Q 2 = 5 GeV 2 . Here, 
F2 = F 2 D + 5F^ + F 2 IP and 8F2 an d F^^, the shadowing corrections to the structure 
functions due to vector mesons and Pomeron-exchange, respectively, are calculated in the 
two phase model. Since the pion component is negligible for Q 2 = 5 GeV 2 , we did not 
include it. 

There is no substantial difference in shadowing between the charm producing (cp) and 
non-charm producing (ncp) parts. (The difference is about 2% in the small x-region). Note 
that the shadowing correction in F 2 p (x, Q 2 ) decreases faster with increasing x, because 
the larger masses of the charmed vector mesons, my, enter in the coherence condition 

2 

T = 'Mxi^ + ( r ^ s ^ ne lifetime of the quark antiquark fluctuation, and M is the 

nucleon mass), compared with the smaller masses of the p and a\. Our results justify the 
assumption that shadowing corrections approximately cancel in the ratio of dimuon and 
single muon cross sections. 



B. Estimate of Parton CSV Contribution 

It would appear that a likely explanation for the deviation of the charge ratio of Eq. |3| 
from one is due to differences between the strange and antistrange quark densities. To test 
this hypothesis, we combined the data in dimuon production, averaged over both neutrino 
and antineutrino events, with the difference between the structure functions in neutrino and 
charged lepton scattering (Eq.([U|)). In combining the neutrino and antineutrino events, one 
measures a flux-weighted average of the strange and antistrange quark distributions. If we 
define a' = N U /(N V + N p ), where N v = 5,030, N D = 1,060 (a' w 0.83) are the number of 
neutrino and antineutrino events of the dimuon production experiment |26|, we have for the 
measured distribution xs(x)^ 

1 1 

xs^(x) = - x [s(x) + s(x)} + -(2a' -l)x [s(x) - s(x)]. (16) 

Now, this equation together with Eq.(]T5]) forms a pair of linear equations which can be 
solved for ^x[s(x) + s(x)] and ^x[s(x) — s(x)]. In this way we can also test the compatibility 
of the two experiments. In addition we have the sum rule that the nucleon contains no net 
strangeness, 

1 

[s(x)-s(x)]dx = (17) 



In the following expressions, we have not enforced the sum rule requirement on the antis- 
trange quark distributions. 

Compatibility of the two experiments requires that physically acceptable solutions for 
^x[s(x) + s(x)] and ~x[s(x) — s(x)}, satisfying both Eq. [15] and Eq. [16], can be found. 
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Using the experimental values a = a' ~ 0.83, we can write x[s(x) — s(x)] = A(x)/8, where 
A(x) = \F^ CFR (x) - 3F^ MC (x) - sw(x), and 5 = (2a - l)/3 « 0.22. Consequently 
even rather small values for A(x) can lead to large differences between s and s. Note that 
the value of the relative neutrino flux, a, depends on the incident neutrino energy. While 
a ~ 0.83 for small x, a is somewhat smaller for higher x-values. However, smaller a would 
lead to an even smaller 5 and would require even larger differences between s and s. 

In Fig. [|we show the results obtained for xs(x) (open circles) and xs(x) (solid circles) 
by solving the resulting linear equations, Eqs. [0| and [16| using the values a = a' = 0.83. 
The results are completely unphysical, since the antistrange quark distribution is negative, 
which is not possible since the distribution is related to a probability. In Fig. [| we show 
the corresponding results for the linear combinations ~x[s(x) + s(x)] (solid circles) and 
^x[s(x) — s(x)] (open circles). The unphysical nature of the solution is demonstrated by the 
fact that \x[s(x) — s(x)] is larger than ~x[s(x) + s(x)]. We also solved the equations using 
the values a = 0.83 and a' = 1 which corresponds to using a subsample of the di-muon data 
containing only neutrino events. In this case, even the sum of the strange and anti strange 
distributions is negative. This is shown in Fig. ^. 

Thus, our analysis strongly suggests that the discrepancy between F 2 CCFR (x, Q 2 ) and 
F 2 (x,Q 2 ) cannot be completely attributed to differences between the strange and an- 
tistrange quark distributions. In other words, assuming parton charge symmetry the two 
experiments are incompatible with each other, even if the antistrange quark distribution is 
allowed to be different from the strange distribution. (Note, that absolutely no restrictions 
were placed on the antistrange quark distribution, aside from the condition that since it 
represents a probability density, it must be non negative.) We stress that our conclusion is 



quite different from that of Brodsky and Ma fl33| , who suggested that allowing s(x) ^ s(x) 



could account for the difference between the two determinations of the strange quark dis- 
tribution. However, they treated the CCFR structure functions as an average beween the 
neutrino and the antineutrino structure function which corresponds to setting a = 0.5. 

At this point there are two possibilities to explain the low-x discrepancy observed between 
the CCFR neutrino and the NMC muon structure functions. Either one of the experimen- 
tal structure functions (or the strange quark distributions) is incorrect at low x, or parton 
charge symmetry is violated in this region, since we have shown that neither neutrino shad- 
owing corrections nor an inequality between strange and antistrange quark distributions can 
explain this experimental anomaly. If we include the possibility of parton CSV, then we 
can combine the dimuon data for the strange quark distribution, Eq. [T^, with the relation 
between neutrino and muon structure functions, Eq. [15], to obtain the relation 

5 -F 2 CCFR (x,Q 2 ) - 3F 2 NMC (x,Q 2 ) - xs^{x) = {2a ~ 1)x [s(x) - s(x)) 
o 3 

+ (3 X [Sd(x) + Sd(x)} - (3 X [Su(x) + Su(x)) 

^ [5u v [x) - 6d v [x)\ (18) 

In Eq. [18] we have used the experimental value a = a', and we have defined the valence quark 
CSV terms Sq v (x) = 5q(x) — 5q(x). We have neglected the effects of possible CSV on the 
extraction of six) from the dimuon data, or on the identification of the structure functions 
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from the neutrino data. This will be discussed below. Since the discrepancy between CCFR 
and NMC data lies primarily in the very small x-region, where the valence quark distribution 
is much smaller than the sea quark, the charge symmetry violation should be predominantly 
in the sea quark distributions. If we set Sq v (x) ~ in this region, Eq. fll8l) can be written as 

\f° ofr {x, Q 2 ) - 3F 2 NMC (x, Q 2 ) - xs^(x) « ( 2a ~ 1)X [s(x) - s(x)} 

+ | [Sd(x) - 5u(x)\ - ^ [Sd(x) + 6u(x)\ . (19) 

Since f3 ~ 0.06 is quite small, CSV arising from the non-isoscalar nature of the iron target 
can be neglected, so in the following we neglect the last term of Eq. [TJ]. 

Using the experimental data we find that the left hand side of Eq. [19| is positive. Conse- 
quently, the smallest value for charge symmetry violation will be obtained if we set s(x) = 
|34| . In Fig.|7| we show the magnitude of charge symmetry violation needed to satisfy the 



experimental values in Eq. The open circles are obtained if we set s(x) = 0, and the 
solid circles result from setting s(x) = s(x). If we use only the neutrino induced di-muon 
events, (i.e we set a' = 1), the coefficient of x[s(x) — s(x)], (5a — 3a' — l)/3, is still positive 
but smaller in magnitude. Consequently, the influence of the uncertainty in s(x) on the 
extracted CSV is smaller. This is shown as open triangles in Fig. |7]. 

In obtaining these results, both the structure functions and the strange quark distribution 
have been integrated over the overlapping kinematic regions and we used the CTEQ4L 
parametrization for f3"5| . In Fig.|S] we show the sensitivity of extracted CSV to the 
parametrization used for s^. The uncertainty due to different parametrizations has been 
partly taken into account since the calculated errors already include the uncertainty of the 
dimuon measurement and most of the parametrizations lie within the experimental errors of 
the dimuon data (except for LO-CCFR s(x)). We note that the magnitude of the observed 
charge symmetry violation in the sea quark distributions is independent of whether we use 
a pure neutrino or antineutrino structure function or a linear combination of neutrino and 
antineutrino structure functions. This is quite different from strange and antistrange quark 
effects which are sensitive to the relative weighting of neutrino and antineutrino events in the 
data sample. Thus, effects due to CSV are independent of the precise value of the relative 
neutrino flux, a. 

The CSV effect required to account for the NMC-CCFR discrepancy is extraordinarily 
large. It is roughly the same size as the strange quark distribution at small x (compare 
the open circles in Fig. [7] with the solid line in Fig. ^). The charge symmetry violation 
necessary to provide agreement with the experimental data is about 25% of the light sea 
quark distributions for x < 0.1. The level of CSV required is two to three orders of magnitude 



larger than the theoretical estimates of charge symmetry violation |36|-|39|1 . Note that, if 



s(x) < s(x) in this region, as suggested in Ref. |33[, we would need an even larger CSV to 
account for the CCFR-NMC discrepancy. 



Theoretical considerations suggest that 5d(x) ~ — 5u(x) |36],|37]]. In fact, since charge 
symmetry violation seems to be surprisingly large, it is reasonable to assume that these 
distributions have opposite signs. We note that with this sign CSV effects also require 
large flavor symmetry violation. One might ask whether such large CSV effects would be 
seen in other experiments. For example, CSV in the nucleon sea could contribute to the 
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observed violation of the Gottfried sum rule ||[37],[3(J and could explain the Fermilab Drell 



Yan experiment ||. This will be discussed in section IV. 

Clearly, CSV effects of this magnitude need further experimental verification. The 
NuTeV-experiment at Fermilab is able to operate either with pure neutrino or pure 



antineutrino beams. The extracted structure functions can be used to build different linear 
combinations, proportional to various combinations of the Sq's and s-s. This will be useful 
to separate CSV from s-s effects. 

At small x, our results can be summarized by 

cT , , ._, s 1/ / s / \ \ 1 fu(x) + d(x)\ 

bd(x) — ou{x) « + s[x)) m - 1 1 



2 v v , v „ 4 y 2 j 

5d{x) + 6u{x) « 0. (20) 

From Eq. |2| we note that such a CSV effect would have little or no effect on the the F2 
structure functions of isoscalar targets, for either neutrinos or antineutrinos. The major 
effect for isoscalar targets would be a significant positive contribution to F$ °(x, Q 2 ) at 
small x, and an equally large negative contribution to F£ °(x, Q 2 ). 

However, if CSV effects of this magnitude are really present at small x, then we should 
include charge symmetry violating amplitudes in parton phenomenology from the outset, 
and re-analyze the extraction of all parton distributions. Given the experimental values 
k = 2S/(U + D) sa 0.5, where S, U and D are the probabilities for strange, up and down 
quarks averaged over x, and the size of CSV effects suggested by the preceeding analysis, 
we would predict that at small x, d n (x) ~ 1.25 u p (x) and u n (x) ~ 0.75 d p (a;). 



IV. EFFECTS OF PARTON CSV ON OTHER OBSERVABLES 

If there is substantial CSV, it should also effect other observables. In the following we 
review the effects which such large CSV terms might have on three quantities; first, the recent 
search for parton "flavor symmetry violation" [FSV] by the Fermilab Drell- Yan experiment 
E866; second, the extraction of the strange quark distribution; and third, experimental 
determination of the Weinberg angle sin 2 (#vy)- 



A. Flavor Symmetry Violation in the Proton Sea 

The results of the recent Fermilab Drell- Yan experiment and the comparison of the 
proton and neutron structure functions measured by the NMC Collaboration Q indicate 
substantial flavor symmetry violation. However, both experimental observations could be 
attributed to charge symmetry violation, as pointed out by Ma || (see also 0). Furthermore, 
both CSV and FSV could be present, as suggested by our analysis of the CCFR-NMC 
discrepancy. Therefore, it is important to examine the effects of CSV on the interpretation 
of the Fermilab and NMC experiments. 

First, we discuss the Drell- Yan experiment which measures the ratio of the dimuon cross 
sections from proton-deuteron and proton-proton scattering. Since CSV is significant in the 
small x region, it is a reasonable first approximation to keep only the contributions to the 
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Drell-Yan cross sections which come from the annihilation of quarks of the projectile and 
antiquarks of the target [41|. In this approximation, the ratio R = a pD / (2a pp ) is given by 



n pD \\ _|_ ^2 _ Sd^l _i_ Ri M _i_ d2 _ 8u 2 l 



4 M 2 

Here, we introduced the notation qj = q(xj) for the quark distributions (x% is the projectile 
x value and 22 refers to the target), and R\ = 

For large Xi?, which corresponds to large xi, the quantity -Ri is small; if we ignore it, we 
have the approximate result 

(22) 

2ctpp 2 1 m 2 1 v ; 

If charge symmetry is valid and c?2 = «2, then we would have R = 1. The experimental 
values give R > 1 at small x 2 ; from Eq. [22|, this could be satisfied if either d 2 > w 2 or <M 2 was 
large and negative. However, the value of 5d(x) extracted from the existing neutrino and 
muon experiments, as discussed in the preceding section, was large and positive at small 
x. The enhancement is on the order of 25% in the small x region where CSV could be 
important. 

In Fig. § the solid circles show the ratio d(x)/u(x) extracted from the Drell-Yan exper- 
iment if we assume the validity of charge symmetry. The open circles in Fig. |9] show the 
result for d(x)/u(x) if we include the CSV term which was extracted from the CCFR-NMC 
data (this is shown in Fig. |7|). Inclusion of parton charge symmetry violation suggested by 
the CCFR-NMC discrepancy plays an important role in the extraction of the FSV ratio 
d/u, in the region x < 0.1. The flavor symmetry violation in the sea has to be substantially 
larger to overcome the CSV term which goes in the opposite direction. In particular, the 
ratio d(x)/u(x) does not approach 1 for small x values. 

We can invert the extracted ratio to obtain the difference [(d — 5d) — u] 

<'"-«3- fi = (':g^+fo -^ +fil - < 23 > 

As a rough approximation, we could neglect 5d in the sum on the right hand side of Eq. [23| 
and keep it in the difference between d and u on the left hand side. For u + d one could use 
a parametrization. This is exactly the way that d — u has been extracted from the Drell-Yan 
data, so that in fact the extracted quantity corresponds to (d — 5d) — u if CSV is present. 

The difference, d — u, can also be extracted from the difference between the proton 
and neutron structure functions measured by the NMC Collaboration || using muon deep 
inelastic scattering. In this case we have 

\{u v {x) - d v (x)) - ^-(Fg(x) - F%{x)) = (d(x) - u(x)) - l(5d(x) + Sd(x)) - hsu(x) + 6u(x)). 
2 2x 3 6 

(24) 

We can make the approximations 6q(x) ~ 5q(x) and 6d(x) ~ — 5u(x), (the latter may not 
be a good approximation since we have FSV), and obtain 
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l -{u v {x) - d v {x)) - ^(F$(x) - F%){x) » [{d{x) - 5d(x)) - u(x)}. (25) 

Comparing this with Eq. [23| we see that, in a first approximation, the quantities extracted 
from the two experiments are the same even if both CSV and FSV are present. However, if 
CSV is present, the term 5 d has to be subtracted from the measured quantity to obtain the 
difference d — u. 

We inverted Eq. |25| by dividing both sides by d — 5d + u = u(r2 + 1), approximating 
d — 5d 4~ u on the left hand side of Eq. |25| by a parametrization of d + u and solving for 
r 2 = d{x-i) / u{x<i) . The structure functions and the parton distribution are integrated for 
each data point over the same Q 2 regions as in the analysis of the charge ratio. The result 
is shown in Fig. ^| as solid triangles. If we subtract the contribution of CSV from the ratio 
r 2 we obtain the result shown as open triangles in Fig. [| We see that charge symmetry 
violation, as suggested by the CCFR-NMC discrepancy, considerably enhances the FSV 
ratio d/u in the region x < 0.1. 

It is interesting to investigate the influence of CSV on the Gottfried sum rule. If both 
CSV and FSV are present the Gottfried sum rule can be expressed as 

S G = - - ~ f 1 dx[d{x) - u{x)\ + I f 1 dx[A8d{x) + 6u(x)\. (26) 
o o Jo 9 Jo 

Now, if 5d(x) ~ — 5u(x) we have 

S ° = I ~ I Jo 1 dx{[d{x) " 6d{x)] ~ U{x)h (27) 

so that, although the CSV suggested by the CCFR experiment does influence the magnitude 
of the extracted FSV, it does not change the experimental value of the Gottfried sum rule 
since the extracted quantities appear in exactly the same form in the Gottfried sum rule as 
in the Drell-Yan and NMC experiments. 



B. Extraction of Strange Quark Distributions 

The differential cross section for the production of opposite-sign dimuons, for neutrino 
and antineutrino deep inelastic scattering from an isoscalar target, are proportional to the 
quark distributions, the CKM-matrix elements, the fragmentation function D(z) of the 
struck quark and the weighted average of the semi-leptonic branching ratios of the charmed 
hadrons B c 

dG{VN lZ ^^ +X) ~ {HO + d(0 - H0W cd \ 2 + 2s(0\V cs \ 2 }D(z)B c (c - fx). (28) 
d^dydz 

For antineutrino scattering the quark distributions should be replaced by the corresponding 
antiquark distributions. In this equation, £ is the rescaling variable defined by £ = x(l + 

2 

|jf), with m c the mass of the produced charm quark. The CCFR-Collaboration used this 
expression together with the parametrization of the quark distributions extracted from their 
structure function data to determine the strange quark distributions. The strange quark 
component of the quark sea was allowed to have a different magnitude and shape from 
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the non-strange component. These two degrees of freedom were parametrized by two free 
parameters k and a, respectively. Further, they treated B c and the mass of the charm quark, 
m c , as free parameters and performed a \ 2 minimization to find the four free parameters 
by fitting to distributions of the measured number densities. We note first, that, provided 
5u = —Sd (see Eqs.(|2]) and fl20|)), charge symmetry violation does not effect the extraction 



of the non-strange parton distributions from the structure function data for small a;-values. 
For an isoscalar target, these distributions can be determined quite accurately, even if charge 
symmetry is broken in the manner given by Eq. However, in extracting the strange quark 
distribution, charge symmetry violation plays an important role. The distribution extracted 
by the CCFR-Collaboration is not the strange quark distribution, but a linear combination 
of the true strange quark distribution and the term in Eq. fl28D coming from CSV. Hence, the 
distribution measured in the experiment, s CCFR (x), is related to the "true" strange quark 
distribution s(x) by 



8{ X ) = S™ F *(x) + ±ffiSm.r). (20) 



2 



cs | 



IV I 2 — 

Since L cd 2 ~ 0.05, the error one makes is roughly two per cent, if 5u is of the same order 
of magnitude as s(x), as the experimental data suggest. 5u(x) is negative and hence the 
true strange quark distribution should be smaller than that determined by CCFR neglecting 
charge symmetry violation. Note that we have neglected all other contributions of CSV to 
the extraction of any other parton distributions, and we neglect higher order corrections, 
which could be sizable • 



C. Determination of sin 2 ( 



It might appear that the precision measurement of sin 2 (9w) from neutrino deep-inelastic 
scattering, carried out by the CCFR Collaboration jyj], rules out the possibility of large CSV 



in the parton sea distributions. Sather [38|] has previously pointed out that the measurement 
of sin 2 (9w) is sensitive to possible CSV effects in parton distributions. 

If charge symmetry is valid the ratio of the differences of neutrino and antineutrino 
neutral-current and charged- current cross sections is given by the Paschos-Wolfenstein rela- 
tion Po" 



„vN -"No i 

R^ S % = \-™\»w). (30) 

°CC ~ a CC z 

The CCFR Collaboration used this relation to extract sin 2 (#n/) and obtained the value 
sm 2 (6 w ) = 0.2255 ± 0.0018(stat) ± 0.0010(syst) |0| which is in very good agreement with 
the Standard Model prediction of 0.2230 ± 0.0004 based on measured Z, W and top masses 
| 32fl . The precision of this result puts strong constraints on CSV in parton distributions. 



However, since the measurement of sin (9yy) based on the Paschos-Wolfenstein relation is 
only sensitive to CSV in valence quark distributions, the substantial charge symmetry 
violation in sea-quark distributions found in this analysis does not contradict the precision 

V)- 



measurement of sin 2 ( 
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This can be seen as follows. The difference between the neutrino and antineutrino 
charged- current cross sections is proportional to the difference between and F 2 and 
to the sum of xF^ and xF%. We see that these linear combinations of the structure func- 
tions are only sensitive to 5q(x) — 5q(x) i.e. CSV in valence quark distributions (see Eq. 

The neutral-current neutrino cross sections on an iso-scalar target, omitting second gen- 
eration quark contributions, is given by 

da uN ° G 2 

j x ^Q2 = 2^ i a « i u ( x ) + d{x) - 5d(x)}x + a d [u(x) + d(x) - 8u(x)]x + 

+ b u [u(x) + d(x) — 5d(x)]x + ba [u(x) + d(x) — 5u(x)]x } (31) 

Here, we defined aj = l 2 +r 2 {l — y) 2 and bf = l 2 (l — y) 2 + r 2 with f — u,d and the couplings 
of the quarks to the neutral-currents are l u = 1/2 — 2/3 sin 2 (6'iy), r u = — 2/3 sin 2 {6w) and 
Id = —1/2 + 1/3 sin 2 (9\y), r d = 1/3 sin 2 (9 W ), respectively. The antineutrino cross section 
can be obtained by interchanging quarks with antiquarks in Eq. We immediately see 
that the difference between neutrino and antineutrino neutral-current cross sections is only 
sensitive to CSV in valence quark distributions. Thus the large CSV effects in the nucleon 
sea quark distributions, suggested by the CCFR-NMC discrepancy, do not influence the 
measurement of sin 2 (9w) based on the Paschos-Wolfenstein relation. 



V. TEST OF PARTON CSV FROM W PRODUCTION AT HADRON COLLIDERS 

Clearly, it is important that the charge symmetry violating distributions, 5d and Su, 
enter with different weights in any observable. Otherwise effects due to CSV are not mea- 
surable. In this connection we also note that most of the measured physical observables are 
proportional to the sum rather than the difference of the charge symmetry violating quark 
distributions. However, 5d and 5u are weighted with the charges of the quarks in electro- 
magnetic interactions, such as deep inelastic scattering with charged leptons and Drell-Yan 
processes. In fact, a comparison between charged lepton and neutrino induced structure 
functions was necessary to detect CSV as we have shown in this paper. We also discussed 
the implications of CSV on the Drell-Yan process. In the following, we show that H^-boson 
production in proton deuteron collisions can also be used to test the CSV found in this 
paper, if we define a suitable observable. Such measurements could be carried out at RHIC 
and LHC. Vigdor originally suggested that asymmetry in H^-boson production could 
be used as a test of parton charge symmetry. 

The cross sections for pD — > W + X and pD — > W~X are given by 

W + X) ~ {m(xi)[-u(x 2 ) + d(x 2 ) — 5m(x 2 )] + d(xi)[u(x 2 ) + d(x 2 ) — 5<i(a;2)]}cos 2 6 c 
+ {u(xi)s(x 2 ) + s(xi)[u(x 2 ) + d(x 2 ) — <M(:r 2 )]}sin 2 c (32) 
W~X) ~ {d(xi)[u(x 2 ) + d(x 2 ) — 5d(x 2 )] + u(xi)[u(x 2 ) + d(x 2 ) — Su(x 2 )}} cos 2 Q c 
+ {u(xi)s(x 2 ) + s(xi)[u(x 2 ) + d(x 2 ) — 5d(x 2 )]}sm 2 O c . (33) 
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We note that the Cabibbo favored terms in the sum of the W + and W~ cross sections are 
invariant under the interchange of X\ and x 2 , if charge symmetry is valid. However, the 
Cabibbo suppressed part of the sum contains terms which are not invariant under X\ X2, 
even if charge symmetry is a good symmetry. Thus, if we define the forward-backward 
asymmetry as 

[XF) (^) w+ M+(^r-(x F )+(^r + (-x F )+(^r-(-x F y 1 } 



we see that it will be proportional to charge symmetry violating terms and terms containing 
strange quarks. Assuming s(x) = s(x) the numerator of Eq. A(-^)(xp) is given by 

A(- — )(xf) = { — [u(x\)Su(x 2 ) + d{x\)5d{x 2 ) + u{x\)8u(x 2 ) + d{x\)Sd{x 2 )\ cos 2 c 
dx p 

+ [s(xi)[d(x 2 ) + d(x 2 ) - Sd(x 2 ) - Sd(x 2 )) sin 2 6 c } - (x\ <-> x 2 ) . (35) 

In the following, we use Su « — 5d, as suggested by our analysis, and note that as the 
charge symmetry violating distribution is of the same order of magnitude as the the strange 
quark distribution, terms proportional to sin 2 c can be neglected. Further, we make the 
approximations Sq(x 2 ) ~ Sq(x 2 ) for x 2 < 0.1 and Sq(xi) ~ for large x\. We then obtain 



Ai 



da 

dxp 



)(x F ) = {[u(xi) + u(xi) — d(xi) — d(xi)]8d(x 2 ) 

+ [8u(xi)u(x 2 ) + Sd(xi)d( y x 2 )]}cos 2 Q c . (36) 

For large xp, the forward-backward asymmetry (due to the first term in Eq. |36|) is pro- 
portional to Sd times the difference between the up and down valence quark distributions. 
The second term is sensitive to CSV in valence quark distributions. However, if Sd ~ —5u 
for valence quarks, as suggested by theoretical considerations [[37]], the second term of Eq. 
|5B] is approximately [d(x 2 ) — u(x 2 )}5d(xi) and is only non-zero if we have FSV. Further, if 



5d{x\) is positive for large x%, as theoretical calculations suggest ||37jj38| then the second 
term will contribute positively to the asymmetry, since d — u > 0, so that it would enhance 
any asymmetry expected on the basis of CSV in the sea quark distributions suggested by 
the NMC-CCFR data. 

We calculated the expected asymmetry A(x F ) for y/s = 500 GeV and y/s = 1000 GeV 
using the values of 5d extracted in section II. The results are shown in Fig. 0. The 
error bars represent the errors associated with Sd and do not include the errors of the W 
experiment. In the calculation, we retained all terms in Eq. [33|. The result obtained by 
using the approximation in Eq. |36] differs only by a few percent from the full calculation. 
We predict considerable asymmetries for large xp. 



VI. CONCLUSIONS 

In conclusion, we have examined in detail the discrepancy at small x between the CCFR 
neutrino and NMC muon structure functions. Assuming that both the structure functions 
and strange quark distributions have been accurately determined in this region, we explored 
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the possible reasons for this discrepancy. First, we re-examined the shadowing corrections 
to neutrino deep inelastic scattering and concluded that shadowing cannot account for more 
than half of the difference between the CCFR and NMC structure functions. Next, we com- 
pared two determinations of the strange quark distributions: the "direct" method, obtained 
by measuring opposite sign dimuon production from neutrino and antineutrino reactions, 
and by comparing the CCFR and NMC structure functions. The strange quark distribu- 
tions extracted by these two methods are incompatible with each other, even if we allow the 
antistrange quark distribution to differ from the strange distribution in an unconstrained 
fashion. 

The only way we can make these data compatible is by assuming charge symmetry viola- 
tion in the sea quark distributions. The CSV amplitudes necessary to obtain agreement with 
experiment are extremely large - they are of the same order of magnitude as the strange 
quark distributions, or roughly 25% the size of the nonstrange sea quark distributions at 
small x. Such CSV contributions are surprisingly large: at least two orders of magnitude 
greater than theoretical predictions of charge symmetry violation. We discussed their influ- 
ence on other observables, such as the FSV ratio measured recently in a proton deuteron 
Drell-Yan experiment, on the Gottfried sum rule and on the experimental determination of 
the Weinberg angle sin 2 W . We showed that such large CSV effects could be tested by mea- 
suring asymmetries in W boson production at hadron colliders such as RHIC or LHC. Such 
experiments could detect sea quark CSV effects, if they were really as large as is suggested 
by current experiments. 
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FIG. 1. The "charge ratio" R c of Eq. y vs. x calculated using CCFR [13] data for neutrino and 



NMC 1 12] data for muon structure functions. Open triangles: no heavy target corrections; open 



circles: v data corrected for heavy target effects using corrections from charged lepton scattering; 
solid circles: v shadowing corrections calculated in the "two phase" model. Both statistical and 
systematic errors are shown. 
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FIG. 2. The differential cross sections for neutrino (solid circles) and antineutrino (open circles) 
deep inelastic scattering as a function of the variable £ 2 = (1 — y) 2 for x = 0.03 and Q 2 = 4 GeV 2 . 
The solid and dotted lines are the results with and without the Callan-Gross relation, respectively. 
The statistical errors are estimated using the experimental fluxes of neutrinos and antineutrinos. 
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FIG. 3. Shadowing corrections in the (a) charm and (b) non-charm producing parts of the 
neutrino structure function, as a function of x, for a fixed Q 2 = 5 GeV 2 . The dashed (dash-dotted) 
lines stand for VMD (Pomeron) contributions. The solid lines represent the total shadowing. 
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FIG. 4. The strange quark distribution x s(x) (open circles) and antistrange distribution x s(x) 
(solid circles) extracted from the CCFR and NMC structure functions. The difference between the 
CCFR neutrino and NMC muon structure functions ~ (see Eq. [0]) is shown 

as solid triangles. The strange quark distribution extracted by CCFR in a LO-analysis [31| is 
shown as solid stars, while that from a NLO-analysis [26| is represented by the solid line, with a 
band indicating ±la uncertainty in the distribution. Statistical and systematic errors are added in 
quadrature. 
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FIG. 5. 2 X [ S ( X ) + (solid circles) and 2 X [ S ( X ) ~ s(x)] (open circles) as extracted from 

the CCFR and NMC structure functions and from the dimuon production data. See Fig. |] for 
the definition of the other quantities. The strange quark distribution extracted by the CCFR 
Collaboration is shown as a solid line with a band indicating ±lcr uncertainty. Statistical and 
systematic errors are added in quadrature. 
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FIG. 6. 2 X [ S ( X ) + s i x )] (solid circles) and 2 x [s(x) — s(x)j (open circles) as extracted from the 
CCFR and NMC structure functions and from the dimuon production data using a' = 1. 
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FIG. 7. Charge symmetry violating distributions x{8d(x) — 5u(x))/2 extracted from the CCFR 
and NMC structure function data and the CCFR dimuon production data under the assumption 
that s(x) = s(x) (solid circles) and s(x) « (open circles) for a' = 0.83, and s(x) = s(x) (solid 
circles) and s(x) ~ (open triangles) for a' = 1. (For the latter only statistical errors are shown.) 
xs(x) at Q 2 = 4 GeV 2 obtained by the CCFR Collaboration in a NLO analysis 26] is shown for 
comparison (solid curve, with la error band). 
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FIG. 8. Uncertainty in the extracted parton CSV term x{5d{x) — 8u{x))/2 due to the 
parametrization used for the dimuon data on charge symmetry violation. Open circles: LO CCFR 
distribution, solid circles: CTEQ4L parton distribution 35]; solid rectangles: CTEQ4D parton 
distribution; solid triangles: NLO CCFR distribution. Here, except for the most "critical" point, 
only statistical errors are shown. 
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FIG. 9. Solid circles: the ratio d(x)/u(x) vs. x, extracted from the Drell-Yan data of FNAL ex- 
periment E866 assuming the validity of charge symmetry. If CS is violated this ratio corresponds 
to (d(x) — 5d(x))/u(x). The result obtained by correcting for CSV is shown as open circles. The 
ratio d{x)/u{x) extracted from the difference of proton and deuteron structure functions measured 
by the NMC group [|| is shown as solid and open triangles, without and with CSV, respectively. 
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FIG. 10. The forward-backward asymmetry for W production, as defined in Eq. [}4|. The solid 
and open triangles are calculated for yfs = 500 GeV and yfs = 1000 GeV, respectively. For Sd, 
the values extracted from the comparison of the NMC and CCFR structure function are used. The 
errors are the errors of 5d and do not include the errors of the W experiment. 
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